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Abstract 

GPS  timing  plays  a  critical  role  in  modern  practice  of  time  errors  estimate  and  synchronization . 
Big  noise  of  the  GPS-based  measured  data  and  inherent  non-stationarity  of  a  time  error  cause 
major  difficulties  here .  In  spite  of  theoretical  separation  of  the  application  fields  for  the  filters 
(stationary  and  non-stationary  signals),  GPS-based  time  error  processes  require  more  explicit 
practical  answer .  Indeed,  what  process  may  be  practically  treated  as  a  stationary  one  and,  to 
opposite ,  how  to  recognize  a  non-stationary  case?  In  this  report  we  answer  these  questions  by 
numerically  and  show  that  for  the  same  transient  time  the  following  filter  should  be  used  to  get 
the  best  accuracy  for  the  known  initial  fractional  frequency  offset  yQ  (time  error  rate)  of  oscil¬ 
lator,  namely  an  average  smoother  for  \yQ  |  <rx,  the  Wiener  filter  for  rx  <  \y0 1  <  r2  >  (ind  the  Kal¬ 
man  filter  for  r2  <\y{)\>  where  rx  and  r2  are  coordinates  dependent  on  the  required  accuracy .  We 

prove  this  conclusion  by  the  example  of  a  time  error  estimate  of  the  rubidium  standard  based  on 
the  reference  timing  signals  of  the  Motorola  GPS  UT+  Oncore  Timing  receiver . 


INTRODUCTION 

GPS  timing  plays  a  critical  role  in  modern  practice  of  time  errors  “on-line”  estimate  and  synchronization. 
Major  difficulties  here  are  caused  mainly  by  big  noise  of  the  measured  data  provided  by  GPS  receiver  and 
inherent  non-stationarity  of  a  time  error.  Statistically,  once  a  random  signal  exhibits  stationary  nature,  then 
the  optimal  Kolmogorov- Wiener  approach  (Wiener  filter  [1])  is  efficiently  used  and,  conversely  ,  the  Kal- 
man-Bucy  technique  (Kalman  filter  [21)  yields  the  best  estimator  for  the  non-stationary  random  signals.  In 
spite  of  both  approaches  evidentally  covering  all  cases,  test  and  measurement  still  use  an  average 
smoother  [3]  owing  to  its  transparency  and  small  variance,  and  despite  of  an  estimate  bias  caused  inevita¬ 
bly  by  non-stationarity. 

Modern  timekeeping  systems  employ  all  three  filters.  Observing  them  even  in  the  past  Proceeding  of 
PTTI’99,  we  realize  that,  for  instance,  to  detect  the  failure  of  a  single  satellite  clock  three  different  space- 


157 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

NOV  2000 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2000  to  00-00-2000 


4.  TITLE  AND  SUBTITLE 

GPS-Based  Time  Error  Estimates  Provided  by  Smoothing,  Wiener,  and 
Kalman  Filters:  A  Comparative  Study 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

Guanajuato  University, Electronics  Dept., Salamanca,  36730,  Mexico, 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS  (ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 

See  also  ADM001509.  32nd  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting,  28-30  Nov  2000, 
Reston,  VA 


14.  ABSTRACT 

see  report 


15.  SUBJECT  TERMS 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

Same  as 

13 

Report  (SAR) 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


segment  timekeeping  subsystems  are  designed  [4].  Two  subsystems  use  a  direct  average  and  the  three-state 
Kalman  filter  is  implemented  in  the  third  case, yielding  the  most  accurate  estimate  as  compared  to  the  aver¬ 
age  filters.  In  the  synchronization  algorithm  for  the  WAAS  network  they  are  based  on  the  two-state  Kal¬ 
man  filter  [5],  providing  a  control  error  of  50ns.  The  same  two-state  filter  was  used  in  the  new  steering 
strategy  for  the  USNO  Master  Clock[6].  Contrarily,  while  steering  the  cesium-based  primary  clock  of  the 
Geo  Uplink  Subsystem  they  employ  an  aveiage  smoother  [7],  Finally,  to  estimate  a  time  error  of  the  master 
clock  for  the  WAAS  test  transmissions  the  algorithm  was  designed  based  on  the  HR  digital  filter  of  the  2nd 
order  [8],  which  may  be  considered  as  a  non-optimal  Wiener  filter. 

The  result  is  following,  for  rather  the  same  quality  master  clocks  and  time  errors  the  different  estimators 
are  used,  and  it  seems  obvious  that  aiming  to  obtain  the  smallest  estimate  and  synchronization  error,  one 
must  follow  the  rules  to  select  the  filter  in  the  optimal  way.  To  work  the  rules  out  for  the  GPS-based  time 
error  processes,  first,  we  must  realize  what  process  may  be  practically  considered  as  a  stationary  one  and, 
conversely,  how  to  recognize  a  non-stationary  case?  Finally,  what  type  of  the  digital  filters  should  be  used 
to  be  the  most  accurate  in  practice  for  the  certain  transient  time  and  known  rate  of  change  of  a  time  error 
caused  by  the  crystal  oscillator,  rubidium  standard,  cesium  standard,  or  even  hydrogen  maser? 

In  this  report  we  answer  the  questions  in  the  following  way.  We  numerically  study  all  three  filtering  algo¬ 
rithms  based  on  an  average  smoother,  the  Wiener  and  Kalman  filters  for  the  same  common  transient  time 
t„  taken  as  an  average  time  of  a  smoother.  We  then  simulate  the  GPS-based  time  error  random  process  with 
a  constant  initial  fractional  frequency  offset  y0  between  reference  and  local  oscillators  and  study  the  filter¬ 
ing  errors  for  the  proper  y0  .  In  this  way  we  determine  the  ranges  for  y0 ,  in  which  each  filter  exhibits  a 

minimal  either  total  (RMS  variance  plus  mean  bias)  or  a  maximum  error.  We  show  that  for  the  same  t„  the 
following  filter  should  be  used  to  get  the  best  accuracy,  depending  on  y0 : 

•  Average  smoother  is  for  the  range  of  |  y0 |  <  ,• , 

•  Wiener  filter  is  for  r,  <  |y0|  <  r2 ,  and 

•  Three-state  Kalman  filter  is  for  t-2  <  |y()|» 

where  rx  and  r2  are  determined  for  the  total  error  as  ru  and  r2(,  and  for  the  maximal  error  as  rIm  and  r2„„  re¬ 
spectively.  Because  processing  time  influences  the  error  strongly  then  we  study  the  errors  for  the  different 
t„,  finding  out  correspondent  dependencies  r!>2(U)  and  presenting  the  simple  approximation  function 
rt  -  af' 5 ,  where  i  =  \t,2t,\m,2m ,  and  a,  is  a  constant. 

We,  finally,  consider  the  example  of  the  filter  selection  to  get  the  most  accurate  estimate  of  the  time  error 
of  the  rubidium  standard  with  known  y0  employing  the  Motorola  UT+  Oncore  receiver.  First,  through  the 

equality  y()  =  a,t~) 5  we  establish  the  critical  transients  for  the  total  and  maximal  errors,  tini ,  tuv  and  t,l]m , 
7  ,  respectively,  and  expect  that  an  average  smoother  will  give  the  smallest  error  for  tu  <tllU  or 

tlr  <  t,nm,  the  Wiener  filter  should  be  the  most  accurate  for  tlrU  <  tlr  <  tll2l  or  tuXm  <  tu  <  tu2m ,  and,  finally, 
the  Kalman  filter  must  be  the  best  for  tll2l  <  tlr  or  tll2m  <  t„  .  As  a  matter  of  the  fact  we  conclude  that  the 
methodology  holds  true  at  least  for  the  considered  case. 

MATHEMATICAL  MODELS  OF  THE  SIGNALS 

Consider  mathematical  presentation  of  a  timing  signal  of  local  oscillator  and  a  noisy  time  error. 
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An  Oscillator  Timing  Signal 

A  model  of  total  instantaneous  phase  <£>(/)  of  oscillator  is  truly  legalized  in  [9]  as 

<t>(0  =  +  y0)t  +  nDvmmt2  +  cp(0,  (1 ) 

where:  <D0  is  initial  phase  offset, 

yo  is  the  fractional  frequency  offset  from  the  nominal  frequency  vnom  (mainly  due  to  finite  fre¬ 
quency  settability  of  the  clock),  y0[  in  ns/hour]  =  3.6x  _y0[in  parts  of  1  O'12], 

D  is  the  linear  fractional  frequency  drift  rate  (basically  representing  oscillator  aging  effects), 

<p(/)  is  the  random  phase  deviation  component. 

While  subtracting  from  (1)  the  same  type  phase  model  of  a  reference  source  and  then  dividing  the  result  by 
2%vnom  > one  comes  to  the  time  error  error  model 


*(0  =  *0  +(fo  -y0re/)t  +  - 


D-D 


^t2  + 


9(0  -9^(0 

2tcv 


(2) 


In  (2)  one  can  also  take  that  all  the  degradation  sources  yo,ef,  Diefi  and  cp,e/0  of  the  reference  source  are 
negligible  as  compared  to  those  of  the  clock  under  test.  As  a  result,  the  x(t )  model  reduces  to  the  practical 
form  of  a  time  error  of  local  oscillator 


D  .2  ,  9(0 


x(t)  =  x0+y0t  +  —t  + 

2  27TV 


(3) 


A  Noisy  Time  Error 


Basically,  we  measure  a  time  error  x(t)  in  discrete  time,  providing  values  xv  for  discrete  time  points  tv  for 
the  constant  time  interval  A  =  tv-  tv.u  where  v  =  0,  1 ,  2,  . . . .  GPS-based  measurements  add  a  noise  to  a  time 
error,  which  has  a  normal  histogram,  thus,  may  be  modeled  as  a  Gaussian  noise.  Both  a  time  error  and  a 
noise  are  summed  (3)  allowing  presentation  of  a  measured  noisy  time  error  (observation)  and  a  clock  state 
with  respect  to  (3)  in  the  matrix  form  as  follows 

=HvA,v+«ov>  (4) 

=  A  v_,A,v_j  +  nXv ,  (5) 

where  £v  is  a  measured  noisy  time  error  (observation),  Xv  is  3-dimensional  oscillator  (clock)  state  vector 
(time  error,  frequency,  and  aging),  Hv  is  1x3  dimensional  measurement  matrix,  Av  is  3x3  dimensional 
clock  state  transition  matrix,  n0v  and  n>_v  are  jointly  independent  white  noises  with  zero  expectations  and 
covariances  Vv  and  Tv  of  3x3  dimension,  respectively, 

Vv=E{nQvnr0v},  (6) 

Vv=E{nXvnl}.  (7) 


In  discrete  time  the  model  (3)  is  transferred  to  the  form  of 

D  ,  2 

=  xv_,  +  yv_, A ~  A  +nxv> 
where  yv  =  yv_t  +Dv^A  +  nyv;  Dv  =  +  nDv ;  and  nxv,  nyv,  and  nDv 


(8) 


are  correspondent  discrete  noises. 
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Correspondingly,  (8)  allows  writing  all  the  matrixes  for  (4)  and  (5)  as 

,  and  Hv  =  [l  0  o]. 


Xv 

"l 

A 

A2  /2 

K  = 

yv 

’  Av  = 

0 

1 

A 

’  nXv  = 

»y, 

A_ 

_° 

0 

1 

_«Dv_ 

present  (4)  in  a  form  of 

and  describe  the  noise  matrix  (7)  as 


=xv  +/*0v, 


'¥v=E{nxyxv}  =  SDvA 


0  0  0 
0  0  0 
0  0  1 


(9) 

(10) 


(11) 


where  SDv  is  two-side  spectral  density  of  a  continuous  white  noise  of  aging  depending  on  A  [10]  and  ex¬ 
pressed  through  the  time  error  noise  straightforward. 


ESTIMATION  OF  A  TIME  ERROR 

Consider  the  three  statistical  algorithms,  namely  an  average  smoothing,  Wiener’s,  and  Kalman’s. 


Average  Smoother 


Average  smoother  allows  a  non-optimal  estimate  xv  of  a  time  error  x„  (in  a  sense  of  a  minimal  RMS  error) 

based  on  an  observation  4v  (10)-  The  algorithm  does  not  require  any  a  priori  knowledge  about  an  oscillator 
state  model,  time  error,  and  even  an  observation  and  is  straightforward 


& 

iV  i~v~N+\ 


(12) 


where  =  0  if  v  <  0,  N  is  number  of  average  points.  According  to  (12),  the  first  estimate  appears  with 
delay  on  JVpoints  of  the  process,  thus,  a  filter  transient  time  equals  ttr  =  A(Af-l)  . 


Wiener  Filter 


In  discrete  time  domain  the  realizable  Wiener  filter  provides  estimate  through  a  convolution  of  its  impulse 
response  hv  and  an  observation  (10) 


i-v-M  +\ 


03) 


where  =  0  if  v  <  0,  M  determines  the  length  of  hv  that  is  taken  to  be  equal  zero  apart  the  time  interval  0, 


...,  (M-  1)A.  To  get  a  minimal  RMS  error  for  estimate,  first,  in  the  tradition  of  Wiener  define  an  optimal 
unrealizable  response  Hok  =  Sxk/(Sxk  +SnQ),  where  Sxk  is  discrete  power  spectral  density  of  a  time  error 

(oscillator  phase)  that  in  spirit  of  Leeson  [11]  is  taken  here  as  Sxk  ~  afk2,  where  A  is  Fourier  frequency;  a 

is  a  constant;  S}1q  is  constant  power  spectral  density  of  a  white  noise;  k  =  0,  . . K-\  ,  and  K  limits  the  length 
of  the  time  error  sequence  taken  at  the  early  stage  to  estimate  the  spectral  densities  with  enough  accuracy. 
Then  use  a  proper  approximating  filter  with  response  Hk  exp(y(j)/t)  and  come  through  the  inverse  discrete 

Fourier  transform  to  the  optimal  impulse  response  hov  of  a  realizable  filter 
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h„=B0ri\Hte»'},  (14) 

where  B  is  a  constant.  In  our  experiment  we  increase  (reduce)  the  level  of  Sn0 ,  increasing  (reducing)  in 
this  way  a  filter  transient  time. 


Kalman  Filter 


The  three-state  Kalman  filter  is  matched  with  a  clock  model  (3)  [12],  allowing  the  following  algorithm  for 
(4)  and  (5)  to  get  estimates  in  discrete  time 


=  Av_,V,  +Kv(^y -HvAv_,V,)  = 


(15) 


where  a  filter  gain  is  defined  as 

kv  =  rvh'UhvRvh7v+fv)-',  (16) 

where  Rv  =  AV-1R,  ,-lAv-l  +  *Fv  is  a  matrix  of  predicted  errors,  and  the  filtering  errors  are  calculated  as 


Rv  -(I~KvHv)Rv 


(17) 


where  I  is  a  unit  matrix.  Transient  of  the  on-line  operating  Kalman  filter  is  due  to  time  expended  to  get 
estimate  (17).  It  may  be  varied  around  the  optimal  value  by  changing  SDv  in  (11),  so  that  since  SDr  rises 

then  tt,  decreases,  and  vice  versa.  The  total  filtering  error  rises  once  transient  is  a  non-optimal. 


NUMERICAL  STUDIES  OF  THE  FILTER  ERRORS 

All  three  filters,  namely  an  average  smoother  (12),  Wiener’s  (13),  and  Kalman’s  (15), are  examined  here 
for  the  same  time  error  process.  To  show  the  effect,  the  noisy  process  (3)  is  simulated  with  variance  cr  = 
40ns  and  with  both  a  stationary  part  of  a  deterministic  function  (0  <  t  <  25  hours),  in  which  case 
x0  =  yQ  =  D  =  0  in  (3),  and  a  non-stationary  part  (25  <  t  hours)  with  yo  =  -2- 1  O'12  =  -7.2  ns/hour  and  D  =  0. 

Because  we  consider  a  transient  time  as  a  principal  performance  of  a  filter,  then,  to  know  trade-off,  we  ob¬ 
tain  the  same  t„  =  lOhours  for  all  three  filters.  While  providing,  the  transient  of  a  smoother  was  evaluated 
by  its  average  time  t,  =  A(iV-l)  and  that  of  the  Wiener  and  Kalman  filter  finished  at  the  level  of  0  9. 

Figure  1  shows  the  simulated  processes  and  estimates  extracted  by  the  filters.  Figure  2  gives  correspondent 
errors  calculated  as  difference  between  estimated  and  simulated  functions 

sv=xv-xv.  (18) 

Just  as  it  had  been  expected  based  on  the  filter  strategies  for  the  dynamic  range  (25hours  <  t),  the  Kalman 
filter  showed  the  smallest  error,  the  Wiener  filter  was  less  accurate,  and,  the  smoother  stayed  hors-concours 
with  its  biggest  error.  Conversely,  for  the  range  of  a  stationary  noisy  error  ( t  <  25hours),  the  smoother  was 
the  best,  the  Wiener  filter  exhibited  more  big  error,  and  the  Kalman  filter  looked  like  the  worst.  Neverthe¬ 
less,  it  is  obviously  speaks  in  favor  of  the  Kalman  filter  that  its  error  remains  say  rather  the  same  for  the 
both  stationary  and  non-stationary  ranges  (Figure  2). 

Excited  by  the  curiosity  of  the  different  filtering  errors  for  the  stationary  and  non-stationary  processes  with 
the  constant  transient,  we  come  to  another  experiment,  while  simulating  only  a  non-stationary  process  and 
evaluating  (1 8)  for  tlr  =  const  and  various  yo  by  the  total  filtering  error 
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(19) 


and  the  maximal  filtering  error 


(20) 


Before  going  on  to  analyze  the  results,  let  us  study  Figure  3,  which  shows  the  total  (a)  and  maximal  ( b ) 
errors  as  functions  of  y0  for  /„.  =  6hours  for  all  three  filters.  There  an  average  smoother  yields  the  smallest 
total  error  for  y0  <  ru  =  1.57x1  O'13,  the  Wiener  filter  is  the  most  accurate  for  ru  =  1.57x1  O'13  <  y0  <  r2,  = 
4  035xl0"13,  and  the  Kalman  filter  is  for  r2t  =  4.035xl0‘13  <  yo-  The  same  filters  provide  the  smallest 
maximal  error  for  the  ranges  of  yo  <  f'u,  =  7.59xl0"13,  rim  =  7.59xl0'13  <y0  <  r2m  =  9.22xl0'13,  and  r2m  = 
9  22x1  O'13  <yo,  respectively.  We  then  estimate  coordinates  r\t ,  r2t ,  nm  ,  and  rlm  by  changing  i„  ,  and  come 
to  the  correspondent  dependences  (Table  1 ). 


Table  1 

Dependencies  of  the  Coordinates  (Figure  3)  on  the  Filter  Transient 


Transient, 

hours 

Measure 

Total  error  coordinates 

Maximal  error  coordinates 

ru 

r it 

r im 

rim 

1  5 

yQ,  ns/hour 

4.032 

11.376 

20.844 

32.782 

Yo 

1.12x10'12 

3.16x10‘12 

5.79x10  2 

9  106x10  17 

30 

1.804 

4.054 

8.28 

11.376 

5.01x1  O'13 

1.126x10‘12 

2.3x10'12 

3.16x1 0‘12 

45 

0.984 

2.07 

3.996 

6.091 

2.733x1  O'13 

5.748x1  O’13 

l.llxlO"12 

1.692x1  O’12 

60 

0.565 

1.453 

2.733 

3.318 

1 ,57x10‘13 

4.035x10'13 

7.593x1 0'13 

9.21 7x10'13 

75 

0.3636 

1.1506 

1.878 

2.107 

1  OlxlO’13 

3. 196x1  O'13 

5.217x10"13 

5.854x1  O’13 

90 

0.3 

0.8028 

1.1088 

1.5329 

8.36x1 0‘14 

2.23x10  13 

3.08x1  O’13 

4.258x10’13 

105 

0.2578 

0.703 

1.036 

1.272 

7.16x10'14 

1.952x1  O’13 

2.878x1  O'13 

3.533x1 0’1 3 

12  0 

0.2542 

0.4896 

1.1614 

1.6718 

7.06x10‘14 

1.36x1  O'13 

3.226x1  O'13 

4.644x1  O'13 

135 

0.2031 

0.4913 

0.7423 

0.9086 

5.643x1 0"14 

1.3647x1013 

2.062x10"13 

2  524x1 0'1 3 

150 

VXV.  ■ 

0.398 

0.7636 

0.7758 

- 

1.106x10  ;3 

2.121x10"13 

2.155x10"13 

165 

0.0295 

0.3207 

0.583 

0.6548 

3.19x10"14 

8.91x10'14 

1.61 9x1  O'13 

1 ,819x10‘13 

180 

0.2497 

0.5526 

0.6073 

- 

6.937x1014 

1.54x10"13 

1.69x1 0‘ 3 

The  curves  provided  in  this  way  (Figure  4)  were  noisy  because  of  different  length-limited  samples  of  the 
simulated  random  process.  Nevertheless,  it  seems  obvious  that  the  following  approximating  function 
G  x  r  5  is  accurate  enough  to  be  used  in  practical  calculations  for  each  coordinate.  We  approximate  those 


as  follows 


ru  =  2.3  x  t~ 


r2l=6.0xt;, 


M  m 


=  1 1  x f"1 5 , and  rlm  =15 xt; 


-1  5 


(21) 

and  use  (21),  while  considering  the  below-given  example  of  the  filter  selection  for  the  GPS-based  time  er¬ 
ror  process  generated  by  the  rubidium  standard. 
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EXAMPLE:  FILTER  SELECTION  FOR  THE  GPS-BASED 
TIME  ERROR  OF  A  RUBIDIUM  STANDARD 

Measurement  of  the  time  error  of  the  rubidium  standard  had  been  carried  out  based  on  the  Motorola  GPS 
Timing  Receiver  Oncore  UT+  with  average  time  A  =  100s  for  about  30  hours  with  the  initial  error  of 
xQ  =  2.1ns  and  offset  y0  =  -4.7  ns/hour  =  -1.3  x  10  12 .  To  separate  the  ranges  for  each  fdter,  substitute  the 
known  yo  for  each  coordinate  in  (21)  and  come  to  the  following  prediction: 

For  the  total  error  119) 

If  tn  <1.455  hours  then  an  average  smoother  should  be  the  most  accurate 
If  1 .455 hours  <  tu  <  2.1 51  hours  then  the  Wiener  filter  should  be  the  most  accurate 
If  2.151hours  <  t,r  then  the  three-state  Kalman  filter  should  be  the  most  accurate 


For  the  maximal  error  (20) 

If  t  <  4.129  hours  then  an  average  smoother  should  be  the  most  accurate 
If  4.129 hours  <  tn  <  5.07 8 hours  then  the  Wiener  filter  should  be  the  most  accurate 
If  5. 01  ih  ours  <  tu  then  the  three-state  Kalman  filter  should  be  the  most  accurate 


Then  tune  the  filters  step  by  step  for  several  transients  to  satisfy  the  above-determined  conditions  and  esti¬ 
mate  total  (19)  and  maximal  (20)  errors  (Table  2). 


Table  2 

Total  and  Maximal  Errors  of  the  Filters  for  the  different  Transient  Times 


ttr, 

Total  error,  ns 

Maximal  error,  ns 

hours 

Smoother 

Wiener 

Kalman 

Smoother 

Wiener 

Kalman 

1  0 

8.367 

8.210 

10  83 

17.204 

16.626 

36  363 

22 

9.429 

8.225 

9  446 

14.447 

14842 

18  088 

30 

10  447 

8.677 

8,279 

15  536 

14,763 

16  850 

45 

13  363 

10  25 

6.622 

19  021 

14  486 

13.805 

60 

16  942 

12.376 

5.306 

26  099 

16  204 

12.759 

An  analysis  of  Table  2  shows  that  just  as  it  had  been  predicted  the  three-state  Kalman  filter  allows  the 
smallest  both  total  error  for  3.0 hours  <ttl.  and  maximal  error  for  4.5 hours  <tn .  The  Wiener  filter  exhibits 

the  smallest  those  errors  for  2.2hours  =  tu  and  3Mours  =  tlr,  respectively.  An  average  smoother  gives  the 
smallest  both  those  errors  for  t,r  <1.0 hours  and  t„  <2.2hours,  respectively.  In  the  range  of  tlr  <  1  .Ohour 
we  watched  also  for  the  small  error  of  the  Wiener  filter.  This  is  because  of  the  limited  processing  sequence 
available  with  small  average  time.  Thus,  we  have  proved  in  this  way  the  above-given  methodology  gener¬ 
alized  by  Figure  4,  except  the  case  of  tt,  =  1. Ohour,  and,  finally,  to  illustrate  the  real  filtering  process,  we 
bringFigures  5 — 8,  those  show  four  cases  of  a  time  error  estimate  provided  by  all  three  filters. 
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CONCLUDING  REMARKS 

We  have  numerically  examined  in  this  report  the  errors  of  the  three  filtering  approaches,  namely  an  aver¬ 
age  smoothing,  Wiener’s,  and  Kalman’s  once  they  are  employed  to  get  “on-line”  time  error  estimate  in  the 
modern  timekeeping  systems.  As  a  major  result  we  present  Figure  4,  which,  at  first,  practically  answer  the 
question  “What  are  stationary  and  non-stationary  time  error  processes?”  separating  space  for  them  left  and 
right,  respectively,  by  correspondent  curves.  And,  it  is  the  most  important,  Figure  4  allows  selection  of  the 
filter  type  for  the  initial  frequency  offset  j0  (rate  of  change  of  a  time  error)  and  the  filter  transient  4.  Based 
on  this,  once  interesting  of  the  maximal  filtering  error,  we  conclude  that  the  following  filter  seems  to  be  the 
most  accurate  depending  on  yo  and  4  ,  namely 

•  For  jo  =  10""  (crystal)  the  three-state  Kalman  filter  is  the  most  accurate  once  4  >  lhour 

•  Once  jo  =  10'12  (crystal  or  rubidium)  then  an  average  smoother  is  accurate  for  4  <  5hours,  the 
Wiener  filter  is  for  5hours  <  4  <  6hours,  and  the  three-state  Kalman  filter  is  for  6hours  <  4 

•  For  jo  <  10'13  (cesium  and  hydrogen)  an  average  smoother  is  accurate  for  4  <  24hours 


The  results  are  readily  extended  to  the  general  case,  including  aging.  Just  account  the  maximally  possible 
frequency  offset  of  your  oscillator  for  the  measurement  (observation)  and  follow  the  above-given  methodol¬ 
ogy,  for  which  more  satisfactory  justification  we  plan  to  revise  the  results  further  analytically. 
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Figure  1.  Simulated  error,  noisy  observation,  and  estimates  provided  by  the  average  smoothing,  Wiener,  and 
three-state  Kalman  filters  for  the  stationary  and  non-stationary  processes 


Figure  2.  Errors  of  the  average  smoothing,  Wiener,  and  three-state  Kalman  filters  in  the  stationary  and  non- 

stationary  ranges  (see  Fig.l) 
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Time  Error,  ns  ^  Time  Error,  ns 


a)  b) 

Figure  6.  Estimates  for  t,r  =  2.2hours  (a)  and  filtering  errors  (b):  the  Wiener  filter  exhibits  the  smallest  error 
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Time  Error,  ns  Time  Error,  ns 


Figure  7.  Estimates  for  ttr  -  3.0hours  (a)  and  filtering  errors  (b)\  the  Kalman  filter  is  the  most  accurate 


a)  b ) 

Figure  8.  Estimates  for  ttr  -  6.0hours  (a)  and  filtering  errors  (b):  the  Kalman  filter  seems  obviously  like  the 

best  estimator 
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